Vortex pump for dilute Bose-Einstein condensates 
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The formation of vortices by topological phase engineering has been realized experimentally to 
create the first two- and four-quantum vortices in dilute atomic Bose-Einstein condensates [A. E. 
Leanhardt et ai, Phys. Rev. Lett. 89, 190403 (2002)]. We consider a similar system, but in addition 
to the lofFe-Pritchard magnetic trap we employ an additional hexapole field. By controlling cyclically 
the strengths of these magnetic fields, we show that a fixed amount of vorticity can be added to the 
condensate in each cycle. In an adiabatic operation of this vortex pump, the appearance of vortices 
into the condensate is interpreted as the accumulation of a local Berry phase. Our design can be 
used as an experimentally realizable vortex source for possible vortex-based applications of dilute 
Bose-Einstein condensates. 
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Quantized vortices manifest the long-range phase co- 
herence of many-particle quantum systems which are de- 
scribed by a complex-valued order parameter field, and 
their existence and stability is intimately related to the 
superfluid properties of the system [IJ. In a loop encir- 
cling the vortex, the phase of the order parameter under- 
goes an integer multiple n of 27r winding. This number n 
is the circulation quantum number of the vortex. 

The first vortices in trapped gaseous Bose-Einstein 
condensates (BECs) of atoms were created by driving 
dynamically a transition between two different hyperfine 
spin states of the condensed atoms with a rotating laser 
field [3] , as suggested in Ref. 3 . Single-quantum vortices 
have also been nucleated by slicing through the conden- 
sate with a focused laser beam moving faster than the 
critical velocity [1], by colliding condensates separated 
by tailored optical potentials [1| , by stirring the conden- 
sate with a laser beam [6*] , and by rotating the condensate 
with an asymmetric trapping potential [Tj. In harmonic 
traps, multi-quantum vortices are unstable against split- 
ting into single-quantum vortices [1, [13], and hence 
rotation of the condensate with frequencies close to the 
harmonic trapping frequency results in triangular vor- 
tex lattices of single-quantum vortices [ll|, • Recently, 
multi-quantum vortices have been created by a coher- 
ent transfer of angular momentum from a photon of a 
Laguerre-Gaussian laser beam to the condensate p^ . 
Furthermore, giant vortices with circulations up to 60 
quanta have been created from such vortex lattices with 
a tightly focused laser removing atoms from the center 
of the condensate for the giant vortex to appear [l3| • 

In contrast to the above-mentioned dynamical meth- 
ods, vortices can also be created with the so-called topo- 
logical phase engineering technique [H, [M [13, [3 ■ Here, 
the phase winding of the vortex is generated by an az- 
imuthally dependent adiabatic turn of the hyperfine spin 
of the condensate. The first two- and four-quantum vor- 
tices were created in dilute BECs with this method by 
Leanhardt et al. [l^ . They reversed the bias field of the 
loffe-Pritchard trap and verified that the resulting two- 
and four-quantum vortex states were consistent with the 



angular momenta of 2h and AH per particle, respectively. 
After the reversal, the bias magnetic field was brought 
back to its initial configuration and the topological un- 
winding of the vortices was observed. 

In this Letter, we show that by controlling the 
strengths of the bias and the quadrupole fields of the 
loffe-Pritchard trap and an additional hexapole mag- 
netic field, vortex pumping can be implemented, i.e., the 
quantum number of the vortex in the condensate is in- 
creased by a constant amount in each pumping cycle. 
The pumping can be carried out either fully adiabatically 
or partly adiabatically and partly instantaneously. The 
operation of the pump is studied analytically within the 
adiabatic approximation, and computationally using the 
Gross-Pitaevskii (GP) equation for hyperfine spin F = 1 
condensates. In this special case, the reversal of the bias 
field in the presence of the hexapole field creates a four- 
quantum vortex. When the bias field is brought back to 
its initial configuration in the presence of the quadrupole 
field, two quanta of vorticity is lost. Thus the pump cre- 
ates two net quanta of vorticity per cycle. To date, the 
vortex pump introduced here is the only known method 
to produce almost any desired amount of vorticity with- 
out particle loss. Furthermore, the achieved vortex densi- 
ties can possibly be high enough for the experimental re- 
alisation of strongly correlated vortex liquid phases [l^l • 

In the adiabatic approximation for an atom with a pos- 
itive Lande g factor, the spin of the condensate is antipar- 
allel to the local magnetic field B — B± + , where B^ 
is the uniform bias field along the positive z-direction 
and the perpendicular part B± is either the quadrupole 
or the hexapole field given in the polar coordinates by 

Bq{r,<P,t) = S;(i)r[cos(</.)i-sin(0)y], (1) 
Bh{r,<t),t) = B'^{t)r[cos{24>)x - sm{2(f))y], (2) 

respectively. Here, x and y are the unit vectors along 
the positive x- and y-directions, respectively. The above 
fields up to quadratic terms in the coordinates can be re- 
alized experimentally with five coils in total and the field 
strengths B'^{t) and -B^(t) can be controlled by changing 
the currents in the coils 
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Let us assume that the adiabatic approximation is 
vahd and the bias field Bz is initially dominating over the 
perpendicular field such that the spin of the condensate 
is aligned with the z-axis, i.e., is in the weak-field-seeking 
state \F,mz = —F), where is the quantum number 
of the spin along the z-direction. Thus the reversal of 
the bias field introduces a spin rotation about the nor- 
mal vector n of the perpendicular field in the xy-plane 
as 



^ — in^-ri 



(3) 



where the dot product is defined as • n — Txnx + 
J-yUy + TzTLz and Ta is the dimensionless spin operator 
in the direction a. From Eqs. H]) and ^ we obtain the 
normal vectors 



nq[(j)) = sm{(j))x + cos{(f))y , 
'nh{4') = sin(20)x 4- cos(20)y. 



(4) 
(5) 



for the quadrupole and hexapole fields, respectively. 
Thus the reversal of the bias field in the presence of the 
hexapole field steers the initial spin state to 



i?^,W|F,m,) = (-1) 



-^^rn.^\F,-mz), (6) 



changing the vorticity of the system by —Arnz- The above 
result can be obtained by expressing the rotation with the 
Euler angles as i?*i^(7r) = i?i(7r/2 - 2(j))Ri,{TT)Ri(2(j) - 

These results on the reversal of the bias field are 
the basis of topological phase engineering to create vor- 
tices [H m, El. In vortex pumping, however, the 
control sequence is to be cyclic, and hence we have to 
turn the spin back to its original position to end the cy- 
cle. If this is carried out with the same perpendicular 
field as the first reversal, Eq. Q implies that the two 
reversals cancel each other, Rn{—7T)Rfi{TT) — I, and the 
topologically created vortex unwinds as was observed in 
Ref. [ 19i] . In general, only one control parameter cannot 
provide geometric effects in cyclic and adiabatic tempo- 
ral evolution. Furthermore, since the spin rotation axis 
in the bias field reversal does not depend on the abso- 
lute strength of the perpendicular field and the reversal 
of the perpendicular field before a subsequent reversal of 
the bias field results in 7r)i?ri(7r) = i?y^(27r) = / for 

bosons, two different perpendicular fields are required for 
a feasible pumping cycle. Such a cycle is shown in Fig. [T] 
for our pumping scheme. In fact, the paths between the 
initial and the middle point of the cycle can be arbitrary 
with the constraint that either B'^ or vanishes always 
in the adiabatic evolution. 

We begin from the ground state of the system which 
corresponds to 



(7) 



where p{r) is the condensate density and Sij the Kro- 
necker delta By adiabatically ramping up the 



hexapole field, reversing the bias field, and ramping down 
the hexapole field, we have performed the spin rotation 
corresponding to Eq. ([6]), and achieved the 4F-quantum 
vortex state e^'^^'^ p^^'f {r)\F, ruz = F), where the su- 
perscript denotes possible changes in the particle den- 
sity due to the emergence of the vortex. Then, we bring 
the spin state back to its initial position by adiabatically 
ramping up the quadrupole field, reversing again the bias 
field, and ramping down the quadrupole field. This cor- 
responds to the spin rotation i?^^ (—''') which removes 2F 
quanta of vorticity. Thus the total effect of n adiabatic 
pumping cycles on the initial ground state is 



(8) 



We observe that our pump increases the vorticity of the 
system by 2F per cycle. Since the particle density in this 
trapping scenario is independent of the azimuthal an- 
gle (f>, the axial angular momentum corresponding to this 
vorticity is (*„|L^|\p„) = 2nNFh, where N = J p{r)dr 
is the total particle number. 
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FIG. 1: (color online). Possible control cycle in the magnetic 
field parameter space for pumping vortices into a spinor Bose- 
Einstein condensate (a). Panel (b) displays an example of the 
temporal changes of the control fields such that the spin at 
a certain distance from the z-axis is reversed with a constant 
speed. The different parts of the control loop in panel (a) are 
marked in the top of panel (b). 

The topological phase engineering method has been ex- 
plained in terms of accumulation of geometric phases [2l| 
for individual spins of the atoms the condensate is com- 



posed of [17|, [19 
result by Berry 



This interpretation is based on the 
2l| that an adiabatically and cyclically 
turned spin acquires a geometric phase equal to —mp^^ 
where is the solid angle covered by the path of the 
spin with the initial state |F, mp). The result of Eq. ^ 
follows from the fact that the solid angle covered by the 
paths of two different spins separated by an angle d in 
the (/) coordinate is 4(? for the reversal of the bias field 
in the presence of the hexapole field. In the case of the 
vortex pump, the result follows directly from Berry's ar- 
gument, since we need only to consider the path traced 
by a single spin. In our pumping scheme, the solid an- 
gle f2 equals twice the angle between the normal vec- 
tors of the hexapole and the quadrupole fields given in 
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Eqs. (HI) and ([5]). This results in the local Berry phase 
of 20, consistent with pumping two quanta of vorticity 
per cycle. Note that the above interpretation of vortex 
pumping in terms of the Berry phase is not directly re- 
lated to the usual Berry phase accumulated to a macro- 
scopic quantum state of a system during adiabatic control 
parameter cycles. 

The adiabatic turning of the condensate spin by the 
magnetic field requires that the magnetic energy domi- 
nates over the kinetic energy at each spatial point and 
that the field is rotated slowly enough to avoid spin 
flips due to Landau-Zener transitions. Thus to guaran- 
tee adiabaticity, the strength of the external magnetic 
field |S(r)| should be large enough in the condensate 
region. However, this condition is not fulfilled on the 
z-axis when the bias field crosses zero, and hence it is 
desirable to avoid particles from entering this area. This 
can be achieved by introducing an additional optical plug 
potential 



(9) 



where A is the strength of the plug and S its width. The 
plug can be realized for pancake-shaped condensates by 
a focused laser beam as was done, e.g., in Ref. [22|. In 
addition, to improve adiabaticity during the magnetic 
field field reversal, the optical plug tends to stabilize the 
created multi-quantum vortex against splitting. 

The magnetic trapping potential arises due to the spa- 
tial dependence of the perpendicular field. Thus the mag- 
netic trap vanishes in the beginning and at the middle 
point of the pumping cycle, which also challenges the 
adiabaticity of the pump. To prevent the condensate 
from expanding, we employ an additional harmonic op- 
tical potential 



We consider a pancake-shaped condensate for which 
the z dependence of the particle density can be taken 
to have th e Gaussian form /i^"-z)^ where ak — 
■y/ h/{mLUk) is the oscillator length in the direction k. 
This approximation is valid for lOz ^ oJr and implies 
that the z dependence of the GP equation can be fac- 
tored out. We measure the length in the units of a,., time 
in 1/uJr, energy in hcoj-i and magnetic field in fkUr/fJ-B: 
and normalize the order parameter to unity. This re- 
sults in a GP equation which depends only on the di- 
mensionless parameters cq = 2^/2TTN{ao + 2a2)/ (Sciz), 
C2 — 2^/2TTN{a2 — ao)/{3az), A/{hujr), and S/ur, together 
with the time-dependent magnetic field. 
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where m is the atomic mass and ujk the frequency of the 
potential along direction k. This kind of potential can be 
realized with laser fields [l| . 

We have studied in detail the temporal evolution of a 
F — 1 condensate during vortex pumping by modeling it 
with the time-dependent Gross-Pitaevskii equation |23l . 

M 



^ 2m 



(11) 



where 4'(r) is a vector representation of the spinor or- 
der parameter {r\'i'), fis the Bohr magneton, gp = 1/2 
the Lande g factor, and cq — 47rfi,^(ao -t- 2a2)/(3m) and 
C2 = 47rfi^(a2 — ao)/(3m) the coupling constants related 
to s-wave scattering lengths oq and 02 of the atoms 
for different spin channels. We study condensates with 
repulsive atom-atom interactions, for which the atom- 
specific constant Cq is positive. The constant C2 is taken 
to be negative for the condensate to be ferromagnetic. 



FIG. 2: Particle density and the complex phase of the |-F = 
l,mz = —1) (left columns) and \F = l,mz = 1) (right 
columns) components of the order parameter at half integer 
multiples of the pumping period T = 360 /uir- The field of 
view is 12ar x 12ar. 

We choose the parameters used in our numerical com- 
putations according to ®^Rb such that C2/C0 = —0.01. 
Furthermore, we set cq = 250. We note that the atom- 
atom interactions do not play a significant role in the 
operation of the vortex pump and qualitatively similar 
results are obtained, e.g., in the non- interacting case 
and for anti-ferromagnetic interactions. The parame- 
ters of the optical plug defined in Eq. ^ are chosen 
as A = lOhujr and 6 = 2ar. The maximum values of 
the magnetic fields are taken to be B^'^^ = AOfkUr/ fJ-B, 

^/max ^ huJr/inBar), and B^™^'^ = hWr / {fiBar) ■ 

Before beginning the pumping cycle, the ground state 
of the condensate is found by a relaxation method with 
the magnetic field parameters set to values correspond- 
ing to t = in Fig.[T] Figure [2] shows the particle density 
and the phase of the relevant components of the con- 
densate during vortex pumping. The accumulation of 
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FIG. 3: Accumulation of axial angular momentum to the sys- 
tem during (a) adiabatic vortex pumping (T = "iGO/ujr) and 
(b) pumping without any optical potentials (T = lO/tJr), for 
which the fraction of atoms in the magnetic trap decreases 
in time as shown. The inset shows the particle density after 
five full cycles of non-adiabatic pumping with the field of view 
12ar X 12ar. 

two quanta of vorticity per cycle is clearly visible in the 
complex phase of the condensate. As more vortices are 
pumped to the BEC, the density depletion along the z- 
axis grows due to increasing core size of the vortex with 
increasing vorticity. The energetic and dynamic stabil- 
ity of the multi-quantum vortex is not guaranteed if the 
core area is much larger than the size of the optical plug, 
and hence initiation of splitting of the multi-quantum 
vortex into single-quantum vortices is observed for the 
ten-quantum vortex shown in Fig. [21 The observed in- 
stability can be overcome by increasing the strength and 
the width of the plug. 

The number of vortices in each spin component is 
quantized, but the axial angular momentum is a continu- 
ous quantity as shown in Fig.[31[a). The figure illustrates 
the increment of vorticity by four during the field reversal 
in the presence of the hexapole field and the decrement 
of it by two when the field is brought back to its ini- 
tial position with the quadrupole field turned on. The 
deviation of the axial angular momentum from 2nhN af- 
ter n pumping cycles is due to excitations away from the 
instantaneous eigenstate the pump is operated in. 

Although the optical potentials given by Eqs. (O 
and pil)) are experimentally realizable using standard 
techniques, they can render the experiment rather com- 



plicated. However, the vortex pump can also be operated 
non-adiabatically without any optical potentials, since 
the spin components which do not follow the magnetic 
field are not trapped, and hence leave the condensate 
region. In fact, no optical fields were employed in the 
experiments of Ref. [19] with the cost of loosing roughly 
half of the atoms during the reversal of the bias field. 
Furthermore, if the maximum value of the bias field is 
much larger than that of the perpendicular field, an in- 
stantaneous switch between the quadrupole and hexapole 
field results in a negligible atom loss. Thus we suggest 
that the first experimental realization of the vortex pump 
could be in a set-up corresponding to the one used in 
Refs. I, H [SI but with three additional coils to gen- 
erate the hexapole field. Our numerical studies on vor- 
tex pumping with this technique show that in fact the 
total atom loss is considerably suppressed by the pres- 
ence of the pumped multi-quantum vortex along the z- 
axis, causing the particle density to vanish there, see 
Fig. El^b). Since the cycle can be traversed much faster 
than in the fully adiabatic case, there is no time for the 
multi-quantum vortex to split, and hence a clean ten- 
quantum vortex is observed in the inset of Fig. EJ^b) after 
five pumping cycles. 

In conclusion, we have shown how to implement a 
vortex pump for spinor Bose-Einstein condensates using 
experimentally feasible techniques. Adiabatic and non- 
adiabatic operation of the pump was demonstrated com- 
putationally in the case of hyperfine spin = 1 by pump- 
ing vortices up to twelve circulation quanta to the con- 
densate. In contrast to the experimentally created giant 
vortex [14| . our scheme produces a pure multi-quantum 
vortex without accompanying single-quantum vortices. 
This vortex pump is an interesting example of adiabatic 
quantum dynamics for which the control parameters of 
the system are varied cyclically but the system does not 
return to its initial eigenspace. Our vortex pump paves 
the way for vortex-based applications of BECs since it 
can work as a vortex source, similar to a spin polarized 
current source in spintronics. 
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